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ABSTRACT 

The viability of advanced viscoplastic models for nonlinear finite element 
analyses of structural components is investigated. Several uniaxial and a 
multiaxial problem are analyzed using the finite element implementation of 
Freed's viscoplastic model. Good agreement between the experimental and 
calculated uniaxial results validates the finite element implementation and 
gives confidence to apply it to more complex multiaxial problems. A comparison 
of results for a sample structural component (the cowl lip of a hypersonic 
engine inlet) with the earlier elastic, elastic-plastic, and elastic-plastic- 
creep analyses available in the literature shows that the elastic-viscoplastic 
analyses yield more reasonable stress and strain distributions. Finally, the 
versatility of the finite-element-based solution technology presented herein is 
demonstrated by applying it to another viscoplastic model. 

INTRODUCTION 

During the last two decades, a number of stress-strain constitutive models 
(called viscoplastic models) have emerged. The motivation in the development 
of these models has been to incorporate observed high-temperature plasticity 
and creep interactions. This is achieved by treating all the inelastic strain 
(plastic, creep, relaxation, etc.) as a single unified quantity, and thus 
automatically including any interactions among them. 
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Viscoplastic models are intended to provide more realistic descriptions of 
high-temperature, time-dependent, inelastic behavior of materials. The attempt 
to make these models realistic by including as much material science as 
possible renders their mathematical framework complicated. The differential 
equations defining the flow and evolutionary laws of the models are, in 
general, highly nonlinear and mathematically stiff. High-temperature 
structural components (such as those used in aerospace or nuclear industries) 
involve complex geometries as well as thermal and mechanical loadings. Closed- 
form analytical solutions for these problems are intractable. Thus, in order 
to use viscoplastic models in nonlinear structural analyses for component 
design, it becomes mandatory to develop and verify suitable numerical solution 
technologies involving, for instance, the finite element method. 

With these requirements in mind, the advantages of employing realistic 
viscoplastic models for nonlinear finite element analyses of structural 
components are explored by analyzing several uniaxial problems and a multiaxial 
problem. The viscoplastic model developed by Freed 1 , and characterized for 
copper, was employed to describe the time-dependent inelastic behavior of the 
material. This viscoplastic model, in its general multiaxial form, was 
implemented in the finite element program MARC 2 . The implementation was 
exercised first on several uniaxial problems involving isothermal and 
nonisothermal cyclic loadings. The results thus obtained were compared with 
the experimental results from Reference 1. This was done to validate the 
finite element implementation. The validated implementation was then applied 
for stress analysis of a sample structural component. The component, called a 
cowl lip, is part of the leading edge of a hypersonic aircraft engine inlet and 
is subjected to severe thermal loading and gradients. Finally, the versatility 
of the finite-element-based solution technology developed for Freed's model was 
determined by using it for another viscoplastic model developed by Robinson^. 
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For completeness, a brief description of Freed's viscoplastic model is provided 
in the following section. 

FREED'S VISCOPLASTIC MODEL 


Following the general mathematical structure adopted by the majority of 
viscoplastic models, Freed's model also incorporates two internal state 
variables, namely, a back stress and a drag strength. The tensorial valued 
back stress accounts for strain- induced or kinematic hardening, whereas 
isotropic hardening in the model is introduced via scalar valued drag 
strength. A small strain and a small displacement formulation are employed. 
It is assumed that there is no explicit coupling between the static and 
dynamic recovery terms in an evolutionary equation for the internal state. 
This imparts a simpler mathematical framework for the model. 

The total strain rate e. . is written as the sum of elastic 
inelastic (including plasticity, creep, relaxation, etc.) e|f?, and thermal 


et!?, rate components. In symbols 
ij. = + el 1 ? + e[*? (i,j = 1,2,3) 


1 i) 11 ij 

Here e denotes the strain, and a dot over a symbol indicates its derivative 
with respect to time t. 

The material is assumed to be isotropic, and Hooke's law is used to 


( 1 ) 


( 2 ) 


relate the elastic strain rate e?^ to the stress rate o^.; that is, 

■eft ±±- 2 - • 2.. . 

e ij E a ij E a kk S ij 

where E is Young's modulus, v is Poisson's ratio, and Sjj is the Kronecker 

delta. Following Einstein's summation convention, the repeated subscripts in 

equation (2) and elsewhere imply summation. 

The deviatoric stress Sjj has the expression 


S.. = a ■ - - r an $ . . 
ij i j 3 kk ij 


( 3 ) 
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With the back stress denoted by Bjj, the effective stress Ejj is 

defined as 

• = S . . — B . . 
i J i J i J 

The constitutive equations of the model are: 

Flow Law 

E 5 


e in = 0Z ZL1 
U 


Jr 


( 4 ) 


(5) 


where 

I = 1 E E 

j 2 2 ^-ij^ij 

The temperature dependence of the model is mainly contained in the 
thermal-diffusivity, function 0, which is defined as 


0 = 


(- fe) = 


exp 


exp - It" 

( m l 


T > 0.5 T 


m 


M 2^ + l 


T < 0.5 T 


m 


( 6 ) 


(7) 


in which equation Q is the activation energy, k is the Boltzmann constant, 

T is the absolute temperature, and T m is the melting point of the material. 

The function Z, called the Zener-Hollomon parameter, has the following 
expressions : 

AF n ; F < 1 

( 8 ) 

A exp[n(F - 1)] ; F > 1 
where A and n are material constants.. The function F is defined as 
J 


Z = 


'2 

F = — 
r D 


(9) 


Here D denotes the drag strength. 
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Evolutionary Laws 


The differential equations governing the growth of the internal state 
variables (i.e., the back stress Bjj and the drag strength D) are 



and 


D = h 


■T ~ er (°) 


( 11 ) 


where 


I 2 = ez 

In equations (10) and (ll), H, h, and L are inelastic material 
constants. The constant L denotes the limiting value of the back stress at 
kinematic saturation!. The recovery function r is defined through the 
following equations: 


r(G) 


o ; 



R(G) ; D > D o 


( 12 ) 


R(G) 


and 


AG n 1 ; G < 1 
A exp[n(G - 1)]/G ; 


G = 



D 


G > 1 


( 13 ) 


( 14 ) 


where S and D 0 are material constants. 

Further, satisfying the dissipativity condition results in the following 
constraint on the function r 


r > 



( 15 ) 
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where 


B 2 " 2 B ij B ij « 16 > 

This condition must be satisfied for the theory to be thermodynamically 
admissible. 

The highly nonlinear and mathematically stiff nature of the constitutive 
equations requires the development of suitable numerical solution technologies 
(involving, for example, the finite element method) to facilitate the 
application of these models in component design. A brief description of the 
finite-element-based solution technology is provided hereunder. 

FINITE ELEMENT SOLUTIONS 

The finite element solutions using Freed's model for both the uniaxial and 
multiaxial cowl lip problem were obtained by using the general purpose finite 
element program MARC. The MARC program was developed expressly for nonlinear 
structural analyses and is well equipped with sophisticated computational 
algorithms and advanced finite element formulations. 

The viscoplastic model of Freed was implemented in MARC through the user 
subroutine HYPELA. The implementation of Freed's model follows essentially the 
same methodology that was utilized earlier by Arya and Kaufman 4 for Robinson's 
mode 1 . 

The mathematically stiff nature of the constitutive equations poses 
significant difficulty in their time integration. Smart time-integration 
strategies are required to integrate these equations efficiently and 
accurately. In the present work, a self-adaptive time integration strategy 
developed by Arya et al.5 was utilized. The strategy is based on the explicit 
Forward Euler method. An explicit method does not require the assembly or 
inversion of Jacobian matrices and thus is more pertinent for problems with 
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large degrees of freedom. The integration strategy used here has been applied 
to other problems in Reference 6, where it was found to work successfully. 
STRESS ANALYSES 

Several uniaxial problems and a multiaxial cowl lip problem were analyzed 
with this finite element implementation. 

The uniaxial loadings shown in Figure 1 consist of mechanical and thermal 
(both isothermal and nonisothermal ) loadings. The nonisothermal loadings 
included in-phase and out-of-phase loadings. The material constants 
characterizing copper, taken from Reference 1, are listed in Table 1. 

For the stress analysis of the multiaxial cowl lip problem, a finite 
element model constructed by Melis and Gladden 7 was utilized. This model, 
shown in Figure 2, is made up of 3294 solid eight-noded brick elements and has 
4760 nodes. The dimensions of the cowl lip are approximately 6 by 1.5 by 
0.25 in. (15.2 by 3.8 by 0.64 cm). However, for the finite element analysis, 
only the 2-in. (5-cm) central portion of the component was modeled to avoid 
di f f icul t-to-quant i fy end effects. A large number of elements were required 
to handle the severity of the thermal loading and gradients. The severity of 
thermal loading, illustrated in Figure 3, can be judged by noticing that the 
temperature at a critical location of the cowl lip rises from 21 to 758 °C in 
only 0.75 sec (Figure 4). However, by taking advantage of the symmetry along 
the thickness, only half of the model was utilized in numerical computations. 
This resulted in a considerable saving in CPU time. 

RESULTS AND DISCUSSION 
Uniaxial Problems 

Predicted hysteresis loops for isothermal and nonisothermal (both in-phase 
and out-of-phase) loadings are shown in Figure 1 and plotted in Figures 5 to 7. 
The strain range for the isothermal loops is 2 percent and the strain rate is 
0.001/sec. The strain range for nonisothermal loops is 0.8 percent and the 
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strain rate is 1 .5xl0 -8 /sec. To facilitate comparison, the experimental 
hysteresis loops are also plotted in these figures. 

Figure 5 exhibits the calculated and experimental hysteresis loops for 
isothermal loadings at temperatures of 22 and 600 °C. A comparison of these 
figures reveals good agreement between the experimental and correlated 
stabilized hysteresis loops obtained with the finite element implementation. 

The predicted and experimental hysteresis loops for the in-phase and out- 
of-phase loadings are depicted in Figures 6 and 7, respectively. The strain 
range for these loops is 0.8 percent and the strain rate is 1 .5xl0 _8 /sec . By 
comparing these figures, it is found that the calculated loops overpredict the 
stress actually seen in the experiments, especially at the cold end of the 
cycle. In general, however, the agreement between the experimental and 
predicted loops is good. This overprediction of stress at the cold ends is 
inherent in the model, and is also seen in the results of Reference 1. 
Multiaxial Cowl Lip Problem 

The steady-state temperature distribution in the cowl lip, obtained from 
Melis and Gladden?, is depicted in Figure 7. The highest temperatures are 
found to occur along the leading edge. A linear interpolation technique was 
employed to estimate the transient temperature distribution in the cowl lip. 
The temperature values thus obtained were used to perform time-dependent 
inelastic stress analyses. . 

To assess the advantages offered by more realistic elastic-viscoplastic 
analysis, it is appropriate to reproduce some earlier results from elastic, 
elastic-plastic, and elastic-plastic-creep analyses of the cowl lip. These 
results, taken from Arya et al. 8 , are shown in Figures 8 to 10. A primary 
creep law was used by Arya et al. for the elastic-plastic-creep analysis. 
Figures 8 to 10 exhibit the various stress distributions during steady-state 
thermal loading across a cross section of the component. 
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Figure 8 presents the elastic stress distribution in the cowl lip. The 
stress values (in magnitude) are seen to be well above the yield strength of 
the copper cowl lip material and are thus unrealistic. The high temperatures 
and thermal gradients along the edge of the cowl lip are expected to cause 
significant inelastic deformation. This consideration led Arya et al.8 to 
perform elastic-plastic and elastic-plastic-creep analyses for the component. 
The results of these analyses are depicted in Figures 9 and 10, which exhibit 
the stress along the edge of the cowl lip. The magnitude of the largest stress 
for both of these analyses is observed to be reduced appreciably from that 
obtained by purely elastic analysis. 

The results of the current elastic-viscoplastic analysis using Freed's 
model are shown in Figures 11 to 13. Figure 11 exhibits the stress 
distribution along the edge at 0.75 sec. The total strain in Figure 12 is the 
sum of elastic, inelastic, and thermal strains. Figure 13 displays the stress 
at 2.25 sec. A comparison of Figures 8 to 11 reveals that the magnitude of 
maximum stress drops continuously from elastic to elastic-plastic to elastic- 
viscoplastic analyses. In Figure 12 the largest total strain is observed to 
occur along the edge of the cowl lip. It is interesting to note that the 
elastic-viscoplastic model is capable of capturing the redistribution of stress 
even during a small period of steady-state thermal loading between 0.75 and 
2.25 sec. This redistribution of stress is apparent on comparing Figure 11 
with Figure 13. The elastic-plastic-creep analysis of Reference 8 was unable 
to depict this redistribution of stress with time. 

ADAPTABILITY OF SOLUTION TECHNOLOGY FOR OTHER VISCOPLASTIC MODELS 

The versatility of the finite-element-based solution technology developed 
during the present work for other viscoplastic models was also investigated. 

It is very encouraging to note that the solution technology easily adapts 
itself for use with other models. As an example, the cowl lip problem was also 
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analyzed by using Robinson's model^. The material constants used in these 
computations are for a copper alloy, Narloy-Z^. 

Some results of this analysis are shown in Figures 14 to 16. Figures 14 
and 15 depict the stress and total strain distributions along the edge of the 
cowl lip at 0.75 sec. The stress along the edge at 2.25 sec is exhibited in 
Figure 16. These figures exhibit trends similar to those observed with Freed's 
model for copper. From Figures 14 and 15, the largest stress and total strain 
are found to occur along the edge of the cowl lip. A comparison of Figure 14 
with Figure 16 clearly shows the capability of this viscoplastic model to also 
predict the redistribution of stress during the steady thermal loading from 
0.75 to 2.25 sec. 

CONCLUSIONS 

A numerical solution technology employing the finite element method was 
developed to perform the nonlinear structural analyses using viscoplastic 
models. The technology is illustrated for Freed's model but is general in 
nature. Several uniaxial problems, as well as a multiaxial cowl lip problem, 
were analyzed. Good agreement between the experimental and predicted 
hysteresis loops for uniaxial thermomechanical loadings confirms the correct 
finite element implementation of the viscoplastic model. By using this 
validated implementation, a finite element analysis is performed for the 
multiaxial cowl lip problem. The advantages of employing realistic 
viscoplastic models are apparent on comparing the results thus obtained with 
the elastic, elastic-plastic, and elastic-plastic-creep analyses. The 
versatility of the finite-element-based solution technology developed herein 
is demonstrated by applying it successfully to another viscoplastic model, 
specifically, that of Robinson. 
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Table 1 Material constants for copper 
[From Reference 1] 


A, s" 1 50 000 000 

D 0 , MPa 1.5 

E, MPa 165 000 - 125 T 

H, MPa 5000 

h, MPa 500 

L, MPa 25 exp(-T/300) 

n 5 

Q, J/mol 200 000 

S, MPa 14.3 

T m , K . 1356 

a, °C -1 . . . 15x10“® + 5xl0“ 9 T 

v 0.34 


Notes : 

1. 1 MPa = 0.145 ksi 

2. Temperature T in degrees 
Celsius 

[°C = (°F - 32 )/l .8) ] 



TIME TIME 


(a) Isothermal loading. 

(b) Nonisothermal (in-phase) loadings. 

(c) Nonisothermal (out-of-phase) loadings. 

Figure 1. — Cyclic thermal and mechanical loading used 
to generate hysteresis loops by the MARC finite element 
code. 
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Figure 2.— Cowl lip finite element model. Elements, 3294; Figure 3.— Thermal loading cycle for cowl lip problem, 

nodes, 4760. 


TEMPERATURE 

°F 13% 125*1 1152 


970 


826 68*1 5*1 1 


399 329 


°C 



758 679 600 521 


M*l1 362 283 20<l 165 



Z X 

Figure 4. — Steady-state temperature distribution in copper cowl lip. Time, 0.75 to 2.25 sec; cross flow, GH 2 coolant. 
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STRESS. KSl 



MECHANICAL STRAIN MECHANICAL STRAIN 


Figure 5. — Predicted stabilized hysteresis loops for copper 
using Freed’s model (isothermal loading). Experimental 
data from Freed and Verrilli (ref. 1). Strain rate, 0.001/sec. 


Figure 6. — In-phase predicted hysteresis loops for copper 
using Freed's model. Experimental data from Freed and 
Verrilli (ref.1). Temperature, 200^ 500 °C; strain rate, 

1 .5x10 _5 /sec. 



-.004 -.002 0 . 002 . 004 

MECHANICAL STRAIN 

Figure 7. — Out-of-phase predicted hysteresis loops for copper 
using Freed's model. Experimental data from Freed and 
Verrilli (ref.1). Temperature, 200** 500 °C; strain rate, 

2.5x1 0~ 5 /sec. 
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STRESS 

KSl 25.0 11.9 -1.32 -14.5 -2 7.7 -40.9 -54.1 -67.2 -80.4 



MPa 172.4 82.0 -9.10 -100.0 -191.0 -282.0 -373.0 -463.0 -554.0 



Figure 8— Stress in Z-direction in copper cowl lip; elastic analysis. Cross flow, GH 2 coolant. 



Figure 9. — Stress in Z-direction in copper cowl lip; elastic-plastic analysis. Cross flow, GH 2 coolant. 


STRESS 

KSI 


MPa 


11.6 9.55 5.48 1.40 -2.67 



80.0 65.8 37.8 9.7 -18.4 


-6.75 -10.8 -14.9 -19.0 



-46.5 -74.5 -102.7 -131.0 


Y 


Z X 

Figure 10. — Stress in Z-direction in copper cowl lip; elastic-plastic-creep analysis. Time, 0.75 sec; cross flow, GH 2 coolant. 
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STRESS 

KSI 9.5 8.2 5.7 3.1 0.5 -2.0 -*.6 -7.2 -9.7 



MPa 65.6 56.7 39.0 21.3 3.65 -14.0 -31.7 -49.4 -67.1 


Figure 1 1 . — Stress in Z-direction in copper cowl lip; elastic-viscoplastic analysis (Freed’s model). Time, 0.75 sec; cross 
flow, GH 2 coolant. 





.00618 .00590 .00534 .00478 .00422 .00366 .00310 .00254 .00198 



Z X 

Figure 12. — Total strain in Z-direction in copper cowl lip; elastic-viscoplastic analysis (Freed's model). Time, 0.75 sec; cross 
flow, GH 2 coolant. 


STRESS 

KSI 9.2 8.0 5.6 3.2 0.8 -1.6 -4.0 -6.4 -8.8 



MPa 63.4 55.1 38.6 22.0 5.37 -11.2 -27.8 -44.4 -61.0 


Y 


Z X 

Figure 13. — Stress in Z-direction in copper oowl lip; elastic-viscoplastic analysis (Freed’s model). Time, 2.25 sec; cross flow, 
GH 2 coolant. 
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STRESS 



Z X 

Figure 1 4-Stress in Z-direction in copper alloy Narloy-Z cowl lip; elastic-viscoplastic analysis (Robinson's model). Time, 
0.75 sec; cross flow, GH 2 coolant. 



Z x 


Figure 15 — Total strain in Z-direction in copper alloy Narloy-Z cowl lip; elastic-viscoplastic analysis (Robinson's model). 
Time, 0.75 sec; cross flow, GH 2 coolant. 


STRESS 


13.7 11.6 7.40 



MPa 34. 5 80.0 51.0 21.9 


-1.04 -5.27 -9.49 -13.7 


-17.9 



-7.2 -36.3 -65.5 -94.5 -123.4 



► 

Z X 


Figure 16— Stress in Z-direction in copper alloy Narloy-Z cowl lip; elastic-viscoplastic analysis (Robinson's model). Time, 
2.25 sec; cross flow, GH 2 coolant. 
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